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At very high energies, the Landau-Pomeranchuk-Migdal effect reduces the cross 
sections for electron bremsstrahlung and photon e+e~ pair production. The frac- 
tional electron energy loss and pair production cross sections drop as the energy 
increases. In contrast, the cross sections for photonuclear interactions grow with en- 
ergy. In solids and liquids, at energies above 10^" eV, photonuclear reactions dom- 
inate, and showers that originate as photons or electrons quickly become hadronic 
showers. These electron-initiated hadronic showers are much shorter (due to the 
absence of the LPM effect), but wider than purely electromagnetic showers would 
be. This change in shape alters the spectrum of the electromagnetic and acoustic 
radiation emitted from the shower. These alterations have important implications 
for existing and planned searches for radiation from i/^ induced showers above 10^'^ 
eV, and some existing limits should be reevaluated. 



1. Introduction 

Although ultra-high energy (UHE) cosmic rays have been studied for many 
years, their origin is still a mystery. Many cosmic-ray models predict signifi- 
cant fluxes of astrophysical neutrinos with energies above 10^° eV. Proposed 
models consider topological defects, superheavy relics of the big bang^ and 
UHE neutrinos as cosmic rays^. Conventional approaches predict that the 
GZK mechanism^ and gamma ray bursts'* produce neutrinos with energies 
above lO^o eV. 

Several groups have searched for radio or acoustic radiation from elec- 
tromagnetic cascades produced by interacting ly^ and have reported upper 
limits on the cosmic flux of i^e (here, i/e includes i^) at energies of 10^" to 
10^^ eV. These searches probe enormous volumes to reach interesting upper 
limits. These limits depend on a good understanding of the cascades that 
are produced in z^e interactions. 
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2. Radio and Acoustic Waves from Showers 

Most of the searches have involved radio waves. The Glue collaboration 
searched for « 2.3 GHz radio waves from the moon using two radio- 
telescopes^. FORTE used satellite-based receivers to search for 30-300 MHz 
radiation from the Greenland ice pack^. Both ANITA^, a balloon-based de- 
tector, and RICE^, a surface antenna array will search for radio waves from 
Vf, cascades in Antartic ice. The SalSA collaboration plans to search from 
radio emission from underground salt domes^. 

Translating these search results into a limit requires a model of the 
electromagnetic radiation produced by the shower. This radiation has been 
evaluated using a 3-d Monte Carlo shower simulations^*^. The calculations 
add the electromagnetic fields from each particle in the shower. This is a 
computationally demanding process which is only practical for relatively 
low energy showers, below 1 PeV. At higher energies, extrapolations are 
usedii'i2. 

The radiation is large when the electromagnetic fields from the different 
particles add coherently and the radiated energy is proportional to the 
square of the shower energy. This happens when the radio wavelength is 
larger than the transverse spread of the shower, seen along the direction 
of propagation. The positive and negative charges cancel but, because 
of positron annihilation, the overall shower contains about 20% more e~ 
than e+. This electron excess produces the coherent radiation. When 
the wavelength is short compared to the lateral size, phase coherence is 
lost. Radiation from the individual particles adds incoherently, producing 
a much smaller signal. The degree of coherence depends on the width 
of the shower; current measurements are largely, but not completely in 
the coherent domain, so the radiation is sensitive to the transverse shower 
spreading. 

The SAUND collaboration has searched for acoustic radiation from 
induced showers, using data from a set of U.S. Navy hydrophones in 

Bermuda-'^^. The acoustic pulse is generated when the electromagnetic cas- 
cade heats the water, causing it to expand rapidly around the cascade. 
The frequencies are low (in the audio range) because the speed of sound 
is so much lower than that of light. The pulse strength is proportional to 
the energy deposition. The acoustic frequencies arc subject to comparable 
coherence conditions as the radio studies. The frequency spectrum of the 
radiation again depends on the shower width. 

None of the calculations used in these studies considers the eff'ect of 
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photonuclcar interactions on electromagnetic showers. Here, we show that 
photonuclear interactions significantly alter Vf. induced showers, and discuss 
how these interactions affect the shower shape and emitted radiation. 



3. Electromagnetic Interactions at High Energies 

In high-energy interactions, the produced electrons receive an average of 
80% of the i/g energy. The remainder is transferred to the target nucleus, 
producing a hadronic shower. 

Electrons with energies E > 100 MeV lose energy largely via 
bremsstrahlung. At somewhat higher energies, electrons lose their en- 
ergy over a distance scale of order Xq, the radiation length. For ice, 
Xo = 36.1g/cm^. The density of ice depends on its composition (mostly air 
content). In Antartica, the ice is covered by a layer of compressed snow. 
For simplicity here, we use a uniform medium with a density of 1 g/cm^ 
(as for water), so Xq = 36.1 cm. Almost all ice is within 10% of this value; 
snow has a somewhat lower density. 

At very high energies, the Landau-Pomeranchuk-Migdal (LPM) effect 
suppresses the bremsstrahhmg of low-energy photons, increasing the dis- 
tance scale. Radiation of photons with energy k from electrons with energy 
E is suppressed when ^^'^^ 

(1) 

where Elpm is a material dependent constant, 

Elpm = « 7.7 TeV/cm • Xq. (2) 



Here m is the mass of the electron. Eg = my/4jT/a = 21.2 MeV, and 
a « 1/137 is the fine structure constant. When E > Elpm, the effective 
radiation length X w X^^J E / Elpm- 

Most shower studies calculate LPM suppression using Migdals 1956 cal- 
culation of this suppression, albeit with some numerical simplifications^^. 
The degree of suppression depends on a variable s 



ELPAlk 



8E{E - fc)^(s) 



where 1 < ^{s) < 2 increases slowly with .s; s — » oo corresponds to no 
suppression, while s — > gives strong LPM suppression. Figure 1 shows 
the differential bremsstrahlung cross sections for different electron energies. 
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Figure 1. The energy weighted differential cross section per radiation length for 
bremsstrahlung in water as a function of y = k/E for electrons with energies of 640 
GeV (top curve), 6.4 TcV, 64 TcV, 640 TcV, 6.4 PcV, 64 PcV and 640 PcV (bottom 
curve). These curves apply for other materials for electron energies of 0.0023, 0.023, 
0.23, 2.3, 23, 230 and 2300 times Ej^pm- 



Figure 2 compares Migdals calculations with data from SLAC experiment 
E-146. The figure shows the photon spectrum from 8 and 25 GeV elec- 
tron beams passing through 3% and 6% Xq aluminum targets^^. When 
k/E < 10"'*, an additional effect, dielectric suppression, further suppresses 
bremsstrahlung^". At the same time transition radiation from the elec- 
tron entering and exiting the target increases the photon flux. An exper- 
iment at CERN has observed the increase in effective radiation length in 
bremsstrahlung from 149 to 287 GeV electrons^^. 

The cross section for pair production may be similarly reduced; when 
the photon energy k is greater than E^pm, the pair production cross section 
is reduced. Figure 3 shows the differential cross section for different photon 
energies. For a given photon energy, symmetric pairs are suppressed the 
most. 

The reduction in bremsstrahlung cross section corresponds to a reduc- 
tion in energy loss (dE/dx) by the electron. Fig. 4 shows the reduction 
in bremsstrahlung dE/dx as a function of incident electron energy. Also 
shown is the reduction in pair production cross section, as a function of 
photon energy. Photons are affected by LPM suppression at much higher 
energies than electrons. 
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Figure 2. Data (points) from 8 and 25 GcV electrons passing througli aluminum targets, 
from the SLAG E-146 collaboration, compared with calculations based on the Betlie- 
Heitler (dashed histogram), and Migdals LPM suppression (dot-dashed histogram), and 
on Migdals calculations with LPM and dielectric suppression (solid histogram). The 
data is binned logarithmically in photon energy, so the Bethe-Heitler spectrum is 
roughly flat. 



4. Photonuclear Interactions 

In contrast to pair production, for fc > 10 TeV, the cross section for pho- 
tonuclear interactions increases with encrgy^^. The dominant contribution 
to the photonuclear cross section (from 'soft' interactions) can be described 
in terms of photon- Pomeron interactions. The cross section for photon- 
proton interactions rises with the 7p center of mass energy W as W'^'^^. 
At very high energies, the photon may also interact directly with a quark 
in the target, gq. Data from HERA on interactions extends up 

to 1^ = 200 GeV^°, equivalent to 20 TeV photons striking stationary pro- 
tons. Direct jq interactions have not clearly been observed at HERA, so 
predictions about this process have significant uncertainties. 

For oxygen, a Glauber calculation accounts for interactions with multi- 
ple nucleons (shadowing). This moderates the ly-dcpendence of the cross 
section. Since oxygen has only 16 nucleons, the number of multiple inter- 
actions is fairly small. Here, we assume that the effect of shadowing in 
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Figure .3. The differential cross section for pair production in water, as a function of 
X = E/k for photons with energies of 64 TeV (top curve), 640 TeV, 6.4 PeV, 64 PeV, 
6400 PcV, 6.4 EeV and 64 EcV (bottom curve). These curves apply for other materials 
for photon energies of 0.23, 2.3, 23, 230, 2300, 23,000 and 230,000 E^pm- 




Electron or Photon Energy/ELpM 



Figure 4. The electron energy loss (dE/dx) for electron bremsstrahlung (solid line) 
and the photon pair conversion cross section (solid line), relative to the Bethe-Heitler 
predictions {i.e. with no LPM suppression) as a function of E/Ei^pm. 

oxygen and direct photon interactions in oxygen and hydrogen cancel each 
other out, so the photon-H20 cross section follows the Pomeron trajectory, 
a « T4^o.i6 rpj^-g jgg^^jg lower total cross section than in the original 
modcl^^, and a slightly higher crossover energy than is given in Rcf. 
For water, this is a reasonable approximation, but, for heavy nuclei, it may 
somewhat underestimate the shadowing and overestimate the cross-section. 

Figure 4 compares these photonuclear cross sections with the pair pro- 
duction cross section in lead and water. In both materials, photonuclear 
interactions dominate for k > 10'^° eV. The cross-over energy is similar for 
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Figure 5. Comparison between the pair production and photonuclcar interaction cross 
sections in (a) lead and (b) water. Above an energy of 10-^'' eV, the photonucleaj: cross 
section is larger than for pair production. 

diverse solids and liquids because the increase in Oypjaf,^ for heavier nuclei 
is cancelled out by the decrease in Elpm as Xq drops. Lunar soil has a 
density of 1.7 g/cm^^, and should have a similar crossover energy. In gasses, 
because of the reduced density, but similar atomic number, the crossover 
point occurs at much higher energies (about 5 x 10^^ eV in air at sea level). 



5. Shower Length 

A fe interaction produces a high-energy electron plus a hadronic shower 
from the struck hadron. Because of the LPM effect, an energetic electron 
will travel a long distance before losing its energy. Figure 4 shows that the 
electron interaction distance can be approximated 



Xe{E)=XoJ^^ (4) 

V i^LPM 

for E > Elpm- Here, we neglect electro- nuclear interactions and direct pair 
production. At sufficiently high energies, these processes will be the domi- 
nant energy loss mechanisms, and electrons will act like muons^^. However, 
these ultra-high energy electron interactions have not yet been studied in 
detail, and so we neglect them here. 

When the LPM effect is strong (i.e. above 10^" eV), the electron trans- 
fers most of its energy to a single photon. For photon energy k ^ Elpm, 
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Figure 6. Diagram of UHE i/g shower development with photonuclear interactions (not 
to scale). 

the photon pair production distance is 



10 times shorter than X,, at the same energy. When the pair production 
cross section falls below the photonuclear cross section, the photon will usu- 
ally interact hadronically. The hadronic interaction length is 1/crp, where 
cr is the photonuclear cross section per molecule (from Fig 4), and p is the 
target density. In water, p — 3.3 x lO'^^ molecules/cm'^, and, at 10^*^ eV, the 
hadronic interaction length in is about 43 m. Figure 6 shows schematically 
how showers evolve with photonuclear interactions. 

We will use a simple model to compare the development of cascades 
with and without photonuclear interactions. Although inaccurate, these 
models are useful for comparison. Purely electromagnetic showers evolve 
via bremsstrahlung and pair production, with each succeeding shower gen- 
eration containing twice as many particles as the last, each with half of the 
energy. Shower development continues until the average particle energy 
drops below the critical energy, Ec- Below this energy, Compton scatter- 
ing becomes more important than bremsstrahlung and the shower quickly 
dissipates its energy into the medium. This occurs after N = ln2(E/Ec) 
generations. In water, Ec = 126 MeV, so for E = 10^° eV, N = 40. With- 
out the LPM effect, a generation occurs in 1 Xq, which is 36.1 cm in water. 
The shower length is 40^0 ~ 15 m. Some authors^°'^^ have used Xoln2 
for the generation length; this leads to showers that are 30% shorter than 
are given here. 

The LPM effect does not affect the number of generations. However, 
the effective radiation length increases rapidly. When the LPM effect is 
large, the effective radiation length for a generation with average particle 
energy Eg is roughly 




(5) 




(6) 
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This length is the geometric average of Xe{E) and X^{k) For each successive 
generation, X decreases by 1/^/2. 

Figure 7 shows the v,. shower length as a function of energy. For purely 
electromagnetic showers, above 10^^ eV, the LPM effect reduces the cross 
sections and the length increases rapidly. This length increase has been 
experimentally observed^*'. The LPM effect also increases the shower-to- 
shower variation^^, complicating measurements. 

Alvarez-Muniz and Zas also studied the length of electromagnetic 
showers"'^^. They defined the length as the distance over which the shower 
has more than a given fraction (10%, 50% or 70%) of the maximum num- 
ber of particles. Above 10^^ eV, where the LPM effect is significant, their 
length scales as E^f'^ . For the 70% fraction, their lengths are slightly larger 
than are given here: 6 m vs. 4.8 m at 1 TeV, and 37 m vs. 20 m at 10^^ 
eV. With a slightly lower containment fraction, the curves would probably 
agree fairly well. 

For hadronic showers, there are no simple parameterizations. Because 
of the higher final state multiplicity (compared to e+e" pairs) and the 
absence of LPM suppression, hadronic showers develop more rapidly than 
electromagnetic showers with similar energies. If ere, we use a simple model 
which generously overestimates the penetration of hadronic showers: we 
treat them as electromagnetic showers, with each generation having only 
twice as many particles as the previous one. Each generation develops over 
a distance A, the hadronic interaction length; A = 83 cm in water. This 
parameterization is shown by the dotted line in Fig. 7. For the energies 
for which data is available^^, this thickness is more than enough for 99% 
containment of hach'onic showers. At high energies, this parameterization 
might overestimate the shower length by 30-70%. Still, at I0^° eV, the 
shower is 3 times shorter than an electromagnetic shower; by f O^'^ eV, the 
difference is a factor of 100. 

In solids, high-energy tt'^ (and some ff) interact before they can de- 
cay. This happens when the decay length, 7/3cr is larger than A. In water, 
7/JcT > A when > 5 PeV. For the r], interactions predominate for E^ > 3 
EeV. So, a 10^" eV hadronic shower develops through several generations 
before a significant electromagnetic component develops. The electromag- 
netic particle energies will be low enough that the LPM effect will be absent. 
Simulations confirm that the LPM effect only occasionally affects hadronic 
showers, when a high energy tt" or 77 decays^^. 

The dashed brown line in Fig. 7 shows the typical length of hybrid 
showers that develop as shown in Fig 6. In F.g 7, electrons take 80% of the 
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Neutrino Energy (eV) 

Figure 7. Shower length as a function of neutrino energy. The soUd blue line is for a 
purely electromagnetic shower, while the dashed brown line is for a hybrid EM/hadronic 
shower. The dotted red line is for a purely hadronic shower. The difference between 
the hybrid and purely hadronic showers is that the hybrid includes the long electron 
path before bremsstrahlung, and the photon path before it interacts hadronically. For 
the hybrid and electromagnetic showers, the electron is assumed to take 80% of the !/e 
energy. 

neutrino energy and photons take 90% of the electron energies. As Fig. 1 
shows, this is a reasonable energy partition. 

Hybrid showers are much longer than purely hadronic showers because 
of the length of the electron and photon tracks. Above an energy of 10^° 
eV, hybrid showers are much shorter than purely electromagnetic showers. 
The presence of photonuclear reactions greatly shortens the Vf^ showers. 

Much of the length comes from the initial electron trajectory. In both 
treatments, electron energy loss due to direct pair production and elec- 
tronuclear interactions is neglect. With a realistic treatment of these ef- 
fects, both electromagnetic and hybrid showers would become shorter, and 
the fractional difference would increase. 

Above 10^'^ eV, the hybrid shower length exceeds 1 km, comparable to 
the typical thickness of the ice or water used for neutrino detection. Even 
if a vertically downward-going neutrino interacts near the target surface, 
the bulk of the signal (which comes from the end of the shower, when the 
number of particles is largest) will be induced in the rock underneath the 
sensitive medium, limiting the sensitivity to near-vertical showers; this may 
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be particularly relevant for SAUND. Above 1 PeV, the earth is opaque to 
neutrinos, so there are no corresponding upward-going neutrinos. 

Some analyses avoid the shower-length problem by considering only ra- 
diation from hadronic cascades from struck nuclei^. Analyses that include 
emission from the 2nd, photon-produced hadronic cascades may find lower 
limits and/or energy thresholds. 



6. Shower Lateral Distributions 

The characteristic lateral spread of electromagnetic shower is given by the 
Molicrc radius, tm = XqEs/Ec = 6.2 cm in watcr^°. For hadronic showers, 
there is no corresponding simple formula for lateral distributions. 

To compare the widths of electromagnetic and hadronic showers in solids 
(liquids should be similar), we consider data from the CERN LAA project^^. 
The collaboration compared the lateral distributions of 5- 150 GcV elec- 
tromagnetic and hadronic showers in a lead/scintillating fiber calorimeter. 
The showers were produced by electron and tt" interactions respectively. 
They modelled the electromagnetic lateral energy density as 

dE _A 1 

~ 7(^.2-^^2)2 ( ) 

where A is the size of the signal, and B is the shower width. They found 
B Pi 2 cm, almost independent of energy. They parameterized the lateral 
distribution of hadronic showers with 2 components: 

dE Bi B2 

— = — e "^/^i + —e 1^^. (8) 
dA r r ^ ' 

where B\ and Bi are the sizes of the two components and Ai and A2 are the 

lateral spreads of the two components. Neither Ai, A2 nor B^jBy varied 

significantly with energy. 

Figure 8 compares these energy depositions (in terms of charge Q de- 
posited in the calorimeter as a function of radius, dQ/dr), for 150 GeV 
electromagnetic and hadronic showers. For electromagnetic showers, half 
of the energy is deposited within a cylinder of 0.9 cm radius; for hadronic 
showers, 2.8 cm is required for the same containment. For 80% containment 
radii of 1.9 cm and 11 cm are required for electromagnetic and hadronic 
showers, respectively. For 90% containment, the radii are 2.8 cm and 21 cm. 
Depending on the containment criteria, the hadronic showers have radii 3 
to 8 times larger than electromagnetic showers. 

The CERN LAA calorimeter showers developed largely in lead; water 
has much lighter nuclei. However, the lateral spread of the shower depends 
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Figure 8. Lateral profiles of charge deposition dQ/dR for 150 GeV electromagnetic 
(dashed line) and hadronic (solid line) showers as observed by the LAA project. 

on the transverse momentum of the particles produced in the interactions; 
the pt distributions should be very similar for lead and water. 

In very high energies hadronic interactions, hard parton-parton interac- 
tions will dominate the cross section. Scattered high transverse momentum 
(Pt) partons will fragment into high px hadrons"^^. At high collision ener- 
gies, these high px hadrons will widen the hadronic showers. We do not 
estimate the magnitude of the increased broadening here, but it could be 
substantial. 

In contrast, the transverse momentum of electromagnetic showers comes 
largely from multiple scattering. The pt from multiple scattering is inde- 
pendent of energy. A small fraction of the pt does come from the pair 
production and bremsstrahlung reactions themselves. When LPM suppres- 
sion is large, the mean opening angle (and, hence, px) in bremsstrahlung 
and pair production increase by a factor S, the suppression factor for the 
angle-independent calculation^^''^^. However, because multiple scattering 
contributes most of the px, the shower width is not significantly affected. 

Because of the hadronic broadening with increasing collision energy, 
the width ratios measured at 150 GeV should be treated as lower bounds 
at 10^° eV. Detailed simulations are needed to develop a better estimate. 
Even though the low-energy part of the shower is largely electromagnetic, 
it retains the lateral spread acquired during it's high-energy hadronic evo- 
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lution; the CERN LAA measurements were sensitive to both the hadronic 
and electromagnetic components of the showers. 

7. Electromagnetic Radiation from Showers 

Coherent radio Cherenkov emission is dominated by the low-energy part of 
the shower, where the net charge excess of over e+ leads to significant 
radiation^^. The frequency spectrum of the radiation depends on the width 
of the shower. An accurate calculation of the frequency spectrum requires 
a full simulation and a detailed calculation. Here, we will consider some 
simple models which should qualitatively illustrate the features which may 
be expected from a full calculation. 

Coherent Cherenkov radio emission occurs when the radio wavelength 
is larger than the shower width; if the wavelength is small compared to 
the width, coherence is lost. Alvarez-Muniz and Zas studied the effect 
of lateral spreading by comparing the radiation from a full 3-dimensional 
electromagnetic cascade calculation of a 10 TeV shower with a simplified 
1-dimensional calculation^^. The importance of the lateral spread depends 
on the observation angle. Coherence is maximal when the observer looks 
along the Cherenkov angle, 9c = SG'^ in ice (and a similar value in the 
lunar regolith^). For frequencies above 500 MHz, Alvarez-Muniz and 
Zas found that the lateral spread significantly reduces the radiation. This 
is within the frequency range explored by most of current experiments; 
these experiments are sensitive to the lateral spreading. 

An observer looking at the Cherenkov angle Oc sees an apparent shower 
thickness given by the lateral spread times 1/ sin(^c') (about 1.21 in water). 
In hybrid showers, the apparent thickness may be several times higher than 
in purely electromagnetic showers; a the lateral spread due to the hadronic 
interactions could have an important effect at frequencies as low as 60-200 
MHz. At higher frequencies, there will be a significant loss of coherence, 
and the radiation will be reduced. 

Buniy and Ralston'^^ found that the radiation from a shower can be 
found using the Fourier transform of the charge distribution. They worked 
in the low frequency limit, where the lateral spread of the form factor had no 
efi'ect. This approximation fails at high frequencies. The region of validity 
is much smaller for hybrid showers than for purely electromagnetic ones. 

The acoustic radiation from a shower comes from what is effectively a 
line source which radiates largely perpendicular to the shower direction; the 
shower looks like an expanding pancake. The acoustic frequency spectrum 



February 2, 2008 15:25 Proceedings Trim Size: 9in x 6in 



showers 



14 

depends on the lateral spread of the shower. For a speed of sound of 
1500 m/s"^^, the maximum frequency for full coherence (i.e. the shower 
is contained within one acoustic wavelength) drops from 75 kHz to 12.5 
kHz when going from electromagnetic to hybrid showers (for 80% lateral 
containment). This covers most of the frequency range studied by SAUND. 

8. Other Implications 

Although photonuclear interactions only dominate above 10^" eV, they may 
affect showers at considerably lower energies, by introducing hadronic com- 
ponents into largely electromagnetic showers. This could affect searches 
for 10^® eV neutrinos produced by cosmic ray interactions with the cosmic 
microwave background^^. 

At intermediate energies (lO""^^ - 10^*^ eV), photonuclear interactions con- 
vert some electromagnetic shower energy into hadrons. For example, a 10^^ 
eV shower evolves through stages that include « 10 10^^ eV particles and 
w 100 10^'^ eV particles, half photons and half electrons. At 10^^ and 10^^ 
eV, the probabilities for photonuclear interactions are about 12% and 3% 
respectively. These probabilities are high enough that most showers with 
energies above 2 x 10^^ eV will have some hadronic component. The pro- 
duced hadrons may decay and introduce a muon content into the shower. 
Although most tt^ and kaons interact before they can decay, charm and 
bottom hadrons produced in the shower may decay semi-leptonically, pro- 
ducing muons. These muon 'tails' could be used to provide fe directional 
information in experiments like IccCubc'^^. 

The initial hadronic shower from the v^. interaction, plus the delayed 
hadronic shower after the electron to photon to hadronic shower conversion 
could mimic a v^- 'double-bang' event^®. The electron plus photon range 
equals 7/3cr for a t (ct = 290 fs) at a neutrino energy around 1 PeV. The 
probability of photonuclear interactions at this energy is small. However, 
there are large fluctuations in r decay length, in energy division (between 
the lepton and the target nucleon) and measurement, and in shower devel- 
opment; all of these may increase the likelihood of misreconstruction and 
consequent misidentification. An accurate estimate of the misidentification 
probability requires detailed simulations. 

9. Uncertainties 

There are significant uncertainties in these calculations. There are ap- 
proximations in Migdals calculations and uncertainties in the additional 
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suppression mechanisms and the photonuclear calculations. 

Migdals calculations do a good job of describing the SLAC E-146 and 
CERN data, with the apparent exception of the E-146 data on carbon 
targets^^. Nevertheless, they have some limitations. Migdal assumed that 
the scattering was Gaussian; A Gaussian distribution considerably imder- 
estimates the number of large-angle Coulomb scatters, and could therefore, 
under-estimate the suppression. Migdal neglected interactions with the 
atomic electrons in the target. This may be important for \ow-Z targets. 

Two newer and more sophisticated calculations, by Zakharov'*" and by 
Baier and Katkov^^, remedy these problems. Both include accurate models 
of Coulomb scattering and account for atomic electrons by using separate 
elastic and inelastic potentials. Both calculations match the experimental 
data. Unfortunately, code for these calculations is not publicly available 
for use in simulations. However, Baier and Katkov give a cross section 
for high-energy (fc ^ Slpm) pair conversion which agrees with Migdals 
calculation to within about 20%, well within the accuracy needed here. 

None of these calculations explicitly consider hydrogen targets. Hydro- 
gen is problematic because the standard Thomas-Fermi screening calcu- 
lations are only accurate for atomic numbers Z > 5. A hydrogen-specific 
screening correction is required to accurately find the cross sections'*^. How- 
ever, because pair production in water is dominated by interactions with 
oxygen, the error in the hydrogen screening causes much less than a 10% 
effect on the cross sections for water. 

Other suppression mechanisms may enter at very high energies. When 
the formation length (reaction zone) is long enough, a nascent photon may 
interact, by either pair production or photonuclear interactions) before it is 
fully formed. Bremsstrahlung and pair production may suppress each other, 
and photonuclear interactions may suppress bremsstrahlung. These effects 
may be important when the formation length (including the LPM effect) is 
larger than Xq or 1/crp. The former can only occur when E > Ep, where 
Ep = 540 TeV in water^^. At this energy, it only applies for a very narrow 
range of k/E (or, for photons, E/k). The range of k/E (or E/k) where it 
applies is only significant at much higher energies: above 10^*^ eV. However, 
this mechanism won't change the conclusion that photonuclear reactions 
dominate above 10^° eV. If anything, the mutual suppression further re- 
duces the electromagnetic cross sections, strengthening this conclusion. 

A broader question for all of the electromagnetic interactions involves 
higher-order reactions like eN e~^e~eN (direct pair production) and 
-yN — > e+e~7^^. Normally, the cross section for these higher order reactions 
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are a factor of order a = 1/137 smaller than the leading-order processes. 
However, the higher-order processes require a higher momentum transfer 
from the target, and so are much less subject to LPM suppression. When 
LPM suppression is large, these higher-order reactions could dominate, cre- 
ating a floor for the electromagnetic cross sections. Detailed comparisons 
have not been done, but 1/137 suppression is only reached at energies above 
1020 eV. 

The final caveat involves photonuclear interactions. The extrapolation 
to 10^" eV is a factor of 2000 in W beyond the HERA data. Between 
the possible moderation of the VF°'^^ Pomeron trajectory and uncertainties 
in direct cross sections, the uncertainties are considerable. However, 
because of the slow energy variation, the large extrapolation uncertainty 
cannot radically change the crossover energy. Even a 50% reduction in the 
rise in cross section above the HERA data (i.e. for A; > 20 TeV) would only 
increase the crossover energy by about 40%. 

In short, although there are significant uncertainties present, even gen- 
erous error estimates do not affect the conclusion that photonuclear cross 
sections are larger than electromagnetic ones for i/g energies above 10^° eV. 



10. Conclusions 

At energies above 10^" eV, photons are more likely to interact hadronically 
than through pair production, and v^. showers are likely to be hadronic. By 
10^^ eV, hadronic showers are 4 times as frequent as electromagnetic show- 
ers. These hadronic showers are considerably shorter, and several times 
wider than purely electromagnetic showers. The increased width of the 
showers reduces the frequencies at which radio and acoustic radiation are 
coherently emitted, and may affect the conclusions of experiments that 
study radio and/or acoustic emission from Vf, showers. 
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